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We investigate the ability to constrain oscillatory features in the primordial power spectrum 
using current and future cosmic microwave background observations. In particular, we study the 
observability of an oscillation arising from imprints of physics at the cutoff energy scale. We perform 
a likelihood analysis on the WMAP data set, and find that the current data set constrains the 
amplitude of the oscillations to be less than 0.77 at 2a, consistent with a power spectrum without 
oscillations. In addition, we investigate the fundamental limitations in the measurement of oscillation 
. . , parameters by studying the constraints from a cosmic variance limited experiment. We find that 

\l ' such an experiment is capable of constraining the amplitude of such oscillations to be below 0.005, 

' implying that reasonable models with cutoff energy scales A > 200Hinfl are unobservable through 

, the microwave background. 

CN , 

^: 

^ ■ I. INTRODUCTION 

cn ; 

' Recent observations of cosmic microwave background (CMB) anisotropics by the Wilkinson Microwave Anisotropy 
, Probe (WMAP) 0, Q have raised the possibility of determining the initial conditions for structure formation with 
^ ' high precision. The most recent data arc consistent with a power-law primordial power spectrum, as predicted by 
most variants of the inflationary paradigm A At the same time, marginal detections of a running of the 

power-law index, as well as glitches in the temperature power spectrum across several multipoles, suggest that the 
' primordial power spectrum may deviate from the power-law form, and the presence of such features may be tested 
with future experiments such as the Planck satellite and CMBPol 

Previous works have taken model independent approaches to the measurement of features in the power spectrum. 
These parametrize the primordial power spectrum in band powers or wavelets |^ ^3 or choose a set of 
window functions to regularize the initial power spectrum |l3l Il4| . An alternate, observation-oriented approach 
involves constructing the principal components of the primordial power spectrum, which are guaranteed to be the 
modes that best constrain deviations from the power-law form [l^- Such general methods are ideal for situations 
O [ where the deviations from the power-law form are unknown. However, if a well-motivated form for the deviations 
^ . exists, the constraints can be improved by incorporating the specific form in the analysis. 

^ ■ A form of modification that is favored by several mechanisms is that of an oscillation on top of the usual power-law 
. . form. Although such a modification can arise from interactions of the inflaton with other fields , we specifically 
. ^ consider the modifications that arise from an ambiguity in the vacuum choice at some high-energy cutoff scale. The 
, presence of a hig h-ei ierg y cutoff in the description of fluctuations during inflation naturally leads to an oscilatory 
''^ ■ modificatio n Jl7l[T^ll9ll20ll2l| . Current data sets allow a very loose limit on a restricted, two-parameter form of the 
' oscillation |22l|. and an oscillation may be favored j23j. although no constraints have been derived on the parameters 
that define the oscillations. In addition, there is some disagreement over whether such effects might be observable by 
future experiments, and parameter forecasts are absent from current literature |22l |24|. Furthermore, the previous 
studies do not account for the full, four-parameter phenomenology that parametrizes our ignorance of the physics at 
the cutoff scale. 

In this work, we quantify the limit on such modifications by using current CMB data to constrain the parameters 
that define the oscillations. Furthermore, we provide robust forecasts for the oscillation parameters by considering 
the constraints from a cosmic variance limited experiment. We briefly review the theoretical framework in Sec. ITTI 
and describe the statistical method in Sec. IIIII We present the bounds from current CMB data in Sec. IIVI and the 
limitations due to cosmic variance in Sec. We conclude in Sec. I VII 
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II. THEORETICAL FRAMEWORK 



In the inflationary scenario, inhomogeneities in the matter density arise due to quantum fluctuations of the inflaton, 
the scalar field that induces inflation. Computations in such a scenario assume that the quantum field theoretic 
description of the inflaton remains valid to arbitrarily small length scales; however, we expect quantum field theory to 
break down at some length scale, so that a complete description of the quantum fluctuations would have to incorporate 
the unknown physics at the cutoff scale. 

In the absence of information concerning the nature of the high-energy cutoff, we assume that the effects of any 
physics beyond the cutoff is equivalent to ambiguities in the choice of vacuum at the cutoff scale. We describe a model 
that parametrizes our ignorance about the cutoff physics, and briefly review the relation between the primordial power 
spectrum from inflation and the cosmic microwave background. 



A. Deviations from scale invariance 



In this section we briefly review the theory of perturbation generation in an inflationary spacetime, and describe a 
generic model for deviations from a scale- invariant power spectrum. A complete treatment can be found in Ref. |2lj |. 
We follow the notation in Ref. except for a change in metric signature. 

We begin by considering a theory of inflation with a minimally coupled scalar field with a potential V{(j)) 



inz. 1 „ 
-^R--{dcl^)^-V{<j^) 
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where nip = (SttG) ^ is the reduced Planck mass. Since the universe is spatially homogenous and isotropic we will 
use the Robertson- Walker metric: 



a^{vi){-dff' + Stjdx^dx^). 



(2) 



Note that the conformal time rj is taken to run from —oo to 0. We assume a spatially flat metric for simplicity. 

If, for a certain range of (p, the potential V^(0) is nearly constant and thus the field is rolling slowly, then the 
scalar field behaves like a cosmological constant and the universe inflates during this phase. In the slow-roll limit, the 
potential can be characterized by the parameters pH 
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with primes denoting derivatives with respect to conformal time. The conditions e, 5 <C 1 imply that the universe is 
inflating, with e = S = corresponding to pure de Sitter expansion. In general, e and S can be functions of time, but 
in the slow-roll approximation we can take them to be constants. 
Perturbing the metric 



g^.^dx^'dx" -> {g^„ + 5g^,^)dx>'dx'' = a'^{r]) [-(1 + 24>)dry2 + (1 - 2<^)S,jdx'dx^] 



and the scalar fleld 
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we compute the action for the perturbations to be |2 

»5'pcrt. — ~ I d X 

Here, z 



(dv) 
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ji^ij^, and d is the ordinary partial derivative over the three spatial indices. The dynamical degree of 

freedom of the perturbations, v = a (^Scj) + (^j^)*^^! is a linear combination of perturbations in the metric and the 
scalar field. 
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Varying the action I?)), we obtain the equation of motion 



-d' \v = Q 

z 



(8) 



for the perturbation. We quantize the field v by promotin g it to an operator v which can be expanded into Fourier 
modes and their associated raising and lowering operators |26l 



(9) 



The raising and lowering operators obey the commutation relation [at, 0^/] = (5(k — k') if the mode functions satisfy 
the following normalization condition 



/ * / * o ■ 



(10) 



Substituting Eq. into the equation of motion we obtain the the equation of motion for the mode function 



Vk = 0. 



(11) 



The quantum field theoretic description is incomplete without a specification for the choice of a vacuum state, 
defined as a state |0) annihilated by the lowering operator (ak|0) — 0). The choice of a vacuum state |0) of the field is 
equivalent to choosing the initial conditions for the mode functions vj. at some time 770 113 ■ In Minkowski spacetime, 
the vacuum choice corresponding to the initial conditions 



Vk 



iVo) = , v'kiijo) = -iVk, 



(12) 



is the unique choice that is invariant under the symmetries of Minkowski spacetime. In curved spacetimes, however, 
no unique choice exists in general. However, the success of quantum field theory in predicting results of accelerator 
experiments implies that length scales much smaller than the curvature of spacetime, the vacuum is essentially the 
same as the Minkowski vacuum. We therefore choose a parametrized set of vacua that reduce to the Minkowski 
vacuum when the curvature H^^ is small compared to the physical length of the mode k. 
Following Ref. [IJ, we impose the initial conditions 
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Here, X and Y are two complex parameters with the normalization condition (|10|l constraining Re(y) to be zero. 
The choice X ^ Y ^ corresponds to the initial conditions for the conventional vacuum choice. 
The parameter 60 is defined as 



aH 



(14) 



which is the ratio of the physical size of a comoving mode k to the Hubble length at time 770. Note that Oq vanishes 
in the limit when the physical size of the mode k, \ — a/k, is much smaller than the Hubble radius (A <C H^^), as 
well as when there is no expansion {H ^0). In these cases, the comoving mode k does not see the curvature of the 
Universe, and we recover the Minkowski vacuum as a result. 

For each mode, there is an ambiguity in the time 770 when Vkirjo) an w^(?7o) £^re specified. Conventionally, the initial 
conditions for all of the modes are specified at r/o —00. However, that limit corresponds to a vacuum choice made 
when the phy sical length of the mode {k/a)~^ — > 0, where we expect conventional physics to break down. Following 
Refs. [T9I I20I l23 | ■ it seems more sensible to define the vacuum for each mode at the time when the physical length of 
the mode equals a fixed length scale A^^, so that 
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This introduces a cutoff scale A into the theory, which can be interpreted as the energy scale where the quantum field 
theoretic description of the perturbations breaks down. The cutoff scale is a free parameter of the theory, although 
this is normally taken to be the Planck scale. With this prescription, the vacuum for each mode k is defined at a 
different time 770(^)7 ^Lud Oq has the following scaling relation: 



X 



(16) 



where iJ, is the Hubble parameter evaluated at the time the physical size of a comoving reference scale fc* crosses the 
cutoff A- 1. 

In Eq. (|13() . X and Y can generally be functions of k. However, since they are dimensionless quantities, they can 
only depend on ratios of dimensionful quantities. Without introducing new scales into the theory, the only such ratios 
are k/aH and higher time derivatives such as {k / a)' / H' . In the former case, we can expand X and Y as functions 
of 6*0 and recover the form of Eqs. In the latter case, we expand X and Y as functions of 6*0, the slow-roll 

parameters, and their time derivatives. Since we have fixed the slow-roll parameters as constants of the theory, they 
can similarly be absorbed into the parameterization without any loss of generality. 

The solution to the equation of motion is given by 



where the Hankel function index is related to the slow-roll parameters by 

= I + 2e - <5. 

The integration constants Ak and Bk are fixed by the initial conditions to be 
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Different choices of vacua leave different imprints upon the temperature anisotropics in the cosmic microwave 
background radiation. For scalar perturbations, the amplitude of the fluctuations can be characterized by the power 
spectrum of the Bardeen variable C = v/ z, defined through the two-point correlation function as 



(0|C(x,ry)C(x + r,77)|0) 



dk sin kr 
k kr 

Substituting Eq. ® into Eq. and using Eq. Q, we find that 

1 fc3 , „ 



A^(fc). 



A (fc) 



(22) 



(23) 



With these equations, we can compute the power spectrum for long wavelength modes \kri\ ^ 1, which is the spectrum 
of fluctuations that seed the microwave background. Substituting Eqs. (|17|l . (|19|l . H2U|) and rt^, = 1 — 4e -I- 25 into Eq. 
H23|l . we obtain the result 



A^(fc) 



^0 



1 + 6Io|A:|cos 



^0(1 -e) 



<y9 — 27r 



Here, is the phase of X, i.e., X = |X|e"^, and Aq is the uncorrected power spectrum 



A^ = 



(24) 



(25) 



with S'q specifying the normalization of the power spectrum. Note that Y does not appear to first order. 
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The modification to the power spectrum (|24|l consists of an oscillatory piece with a fc-dependent amplitude. A more 
convenient choice of parameters would be related directly to the characteristics of the oscillation, rather than to the 
theoretical model. To rewrite the power spectrum, we introduce the following parameters: 



2e A 



3 - ris 



(26) 
(27) 

(28) 



The parameter A is the amplitude of the oscillation, modulo the k dependence (k/k^)~'^, and lo is the oscillation 
frequency in log(fc/fc,). Lastly, a characterizes the phase of the oscillation. Previous studies [22, 113, 01 do not 
explicitly include the phase as a free parameter in their computation of constraints. However, since the phase shifts 
the starting point of the oscillations, this should not be neglected in assessing the observability of the oscillations. 

With these parameters, the power spectrum (|24|) is rewritten in its final form 



Alik) 



1 + A 



cos 



(29) 



The choice X = uj = 2e/C, and a = — 7r/2 corresponds to the two-parameter model considered in '22*1, with 

The amplitude of the oscillation, |A|-^, is dependent on the parameters |A| and A, but the frequency jtij"^ 
depends only on A. Thus a measurement of the frequency u allows us to break the degeneracy in the amplitude. 

Finally we comment on the possible range of values that \X\ and \Y\ can take. Because Eqs. ()13|l are expansions 
in ^0, there is no upper bound for \X\ and \Y\ if we do not require the expansion to converge. However, there are a 
couple of reasons why we should expect |A"| and |y| to be small. The first is a theoretical prejudice: if \X\ and |F| 
are not small, then there is no reason to expect the 0{9q) and hi ghe r-order terms to be negligible. In addition, limits 
on gravitational particle production may constrain |A| < 5.5^ |2j|. 

B. Primordial Power Spectrum to the CMB 

The temperature and polarization anisotropics of the microwave background arise from acoustic oscillations of 
the inhomogeneous photon-baryon fluid, which are seeded by the metric fluctuations created during inflation. The 
multipole moments of the fluctuations can be related to the Bardeen variable ^(k) as 



47r 



d^k 
(2^ 



r;-^(fc)yr*(k)C(k), 



(30) 



where X = {T,E} is the temperature or E-mode multipole moments, and T^{k) denotes the transfer function 
encoding the acoustic physics and the projection effects (see, e.g., Ref. for discussion on the structure of the 
transfer functions). 

The CMB power spectra, defined as 



(31) 



are then related to the primordial power spectrum A?(fc) by 



/(^ + 1) 
2n 



C, 
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T{'{k)Ti^\k)Ai{k) 



(32) 



The transfer functions in Eq. H32|l limit our ability to measure A^(fc) in several ways. First, (k) projects the 
fiuctuations onto a sphere, and this projection peaks along /peak ~ kr^ with a width of Al w 20, where r» = 13.8Gpc 
is the angular diameter distance to the last scattering surface. Since a Fourier mode with wave number k projected 
onto a sphere also contributes to multipole moments with I <^ ^pcaki power at wave number k gets distributed to 
multipole moments with / <C /peak, as well as in the narrow band with Al « 20 around / <C /peak- 
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The transfer functions T^{k) also include contributions that arise due to secondary effects. In the low multipole 
regime, the temperature transfer function includes a contribution from the integrated Sachs- Wolfe effect [l^, which 
provides additional power at lower multipoles. Rescattering of the CMB photons during reionization contributes to 
a feature in the low multipole regime of the polarization transfer function. This shows up as a bump in the EE 
power spectrum and TE cross spectrum, of which the latter was first detected by the WMAP collaboration 30] . Both 
features tends to broaden features in the primordial power spectrum. Although there are other foreground effects at 
high multipoles, we work in the multipole regime below which those secondaries are insignificant. 

An additional limitation arises due to gravitational lensing of the CMB 31, 32, 33). Lensing has the effect of 
redistributing power between adjacent multipole bins, thereby reducing the width of resolvable features to roughly 
A/ ~ 60. This translates to a maximum log frequency of roughly a; ^ 30, and potentially an important limitation on 
measuring the properties of the primordial fluctuations |l5j . We will restrict the range of lo to be below 30 in our 
analysis to avoid the regime where contamination due to lensing is an important effect. 



We are interested in how well the microwave background can in principle distinguish signatures of new physics 
from the standard inflationary models. Previous works have estimated the uncertainties on parameters related to the 
cutoff physics using the Fisher matrix formalism j24 | , or used a likelihood approach on a fiducial model containing a 
modulation p3 |. Both start from a two-parameter model for the primordial power spectrum pjj. and therefore do 
not incorporate the full range of effects that might arise from the existence of a high-energy cutoff. Furthermore, the 
analyses of the ideal case depend on the choice of fiducial models. We currently have no a priori knowledge, either 
from experiments or theoretical prejudice, about what set of parameters would constitute a reasonable fiducial model 
other than the standard inflationary model without oscillations. However, Refs. [23| and [23| both assume fiducial 
models containing oscillations, and the resulting parameter forecasts therefore only apply to those fiducial models. 

In this work, we instead present two analyses, both of which complement the existing literature. In the first, we 
use current CMB data sets to place constraints on the four-parameter description of the oscillation induced by the 
phenomenology at the cutoff scale. In the second, we present the fundamental limitations on measurements of the 
oscillation parameters that arise from the statistics of CMB measurements. 

Both of these analyses rely on the calculation of the likelihood over the expanded parameter set, which is ac- 
complished using Markov Chain Monte Carlo (MCMC) methods. In the following sections, we briefly describe the 
calculation of likelihoods for microwave background data, as well as details concerning our use of MCMC. 



Let us suppose that we are given a data set C*/, along with a set of model power spectra C/(p), with p denoting a 
vector in an TV-dimcnsional parameter space. The posterior probability distribution P(p|C;), which is the probability 
density at point p in the parameter space given the data set, is then computed by 



Here, P(p) represents our prior assumptions about the parameters, and C denotes the likelihood function. The 
proportionality constant is set by requiring that 



over the allowed range in parameter space. In many instances, we do not have any prior knowledge about what the 
parameter should be, so that we assign a flat prior P{p^) = 1 for the zth parameter. In such a case, the posterior 
probability is proportional to the likelihood function, so maximizing the likelihood function in the parameters is 
equivalent to finding the point with the maximum probability density. To find the la uncertainties around the 
maximum likelihood point, we then find contours that enclose 68.3% of the posterior probability. 

Assuming that the CMB multipole moments are Gaussian distributed, the likelihood for an all-sky, noiseless exper- 
iment is given by 



III. ANALYSIS 



A. Likelihood Analysis 



P(p|(70(xF(p)£(Q|p). 



(33) 




(34) 




(35) 



7 



where the data vector is given by d[" = {T™, EJ^}, and the covariance matrix is 
with C^'^ denoting the theoretical power spectra. Defining the estimator for the power spectra according to 



C, 



XX' _ \x*Txr\ 



21 + 1 

in 

the log likelihood function becomes (up to an irrelevant constant) 



(36) 
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(38) 



We have chosen the constants so that — 21n£ = when the theoretical power spectra C^^ are equal to the data 
estimates , such that — 21n£ behaves as an effective chi-squared statistic. We will use this form of the likelihood 

function in our analysis of the idealized case. 

If only the temperature power spectrum is measured, then the log likelihood function is given by 



2 In/: 



TT 



TT 



(39) 



For a realistic experiment, the presence of sky cuts and instrumental noise necessitate the use of various approximations 
to the likehhoods, Eqs. and 

The transfer functions {k) used to calculate the model power spectra C;^"^ can be obtained by using publicly 
available Boltzmann integral codes such as CMBFAST 34] and CAMB JiSj. However, these codes sample the transfer 
function sparsely in k and I to reduce computation time, and the resulting power spectra can deviate from the true 
power spectra if there are features in the initial power spectrum which have scales comparable to the gridding. In our 
parametrization, higher values of uj correspond to higher oscillation frequencies, so we expect the computation from 
fast Boltzmann codes to become less reliable in that regime. 

To illustrate this, we compare the power spectra generated using CAMB to those from a Boltzmann hierarchy code 
developed by Hu |l3 , which has the advantage of capturing the high frequency oscillations that might be problematic 
with the integral codes. To determine the frequency limit imposed by the accuracy of the code, we compare the power 
spectra generated by CAMB to that from the hierarchy code for various values of lo. As can be seen in Fig. ^ where 
we plot the fractional deviation 



AG 

Ci 



^CAMB _ ^truc 
^true ' 



(40) 



the amplitude of the deviation increases with w. We impose the fairly strict condition that deviations be less than 
0.2%. As a result, the range of lo is restricted to be less than 20 for computations performed using CAMB. 

Since computing the likelihood function from T^{k) on a grid over the entire parameter space is generally not 
feasible for more than two parameters, we utilize Markov Chain Monte Carlo sampling to compute the posterior 
probability distribution. 



B. Markov Chain Monte Carlo 



Markov Chain Monte Carlo (MCMC) methods have been applied to CMB analysis of the data from several exper- 
iments [s^ls^l- We give a brief overview of the method, as well as a description of our implementation, and refer the 
reader to Refs. [36ll38l| for full treatments. 

MCMC is a method designed to simulate the posterior probability distribution through a series of random draws. 
In our case, the MCMC produces a set of N points {pt,t = 1 . . .N} that fairly sample the distribution P(p|C';), so 
that the density of points {pt} is proportional to the posterior distribution. Statistical properties of the parameters 
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can then be estimated using the corresponding properties of the samples; for instance, the mean and variance of 
parameter can be estimated by using the sample mean and variance via 

1 ^ 



N 
t=i 



N 



'''^p'^ = wi^T.^p^-p'^'- (41) 

Our goal is to produce a chain whose stationary distribution is the distribution P{p\Ci). We describe the Metropolis- 
Hastings algorithm js^ , used in many MCMC implementations for CMB analysis including the WMAP analy- 
sis and the public package CosmoMC [sBl, as well as in our own implementation. We refer the reader to Ref. 
for details on CosmoMC, which we utilize in Sec. Irvl and give a generic overview of the algorithm here. 

At each step t in this algorithm, a candidate point q is first chosen from a proposal distribution /3(q|pt). The point 
q is accepted with probability 

. ( P(q|Q)/3(p,|q) ^ 
7r(Pt, q) = mm 1, . (42) 

If q is accepted, then the chain continues, with pt+i = q; otherwise, Pt+i — Pt, i.e., the chain does not move. Note 
that the proposal distribution can depend on the last point of the chain. We take /3(q|pt) to be a Gaussian distribution 
with mean pt and a fixed, diagonal covariance matrix. We first run a short Markov chain to determine the variance 
to be used in the proposal distribution. 

Proper conclusions from chains are not possible without determining if the chains have converged to the stationary 
distribution. The first step in ensuring this is to remove the first few thousand points belonging to a burn-in period. 
All points past the burn-in period are assumed to have forgotten about its starting point, and can be thought of as 
being sampled from the stationary distribution. 

Additionallj^ we employ a convergence criterion, advocated by Gelman and Rubin and used by the WMAP 
collaboration |1,|3j|33) to ensure proper mixing and convergence of chains. In this criterion, one generates several 
chains from different starting points, and quantifies whether these chains are indistinguishable from each other. 

Suppose one generates M chains from different starting points, each chain with 2N points. For each parameter p*, 
we take the last N points from each chain, and label the sequences as {p^^}, with a = 1, . . . ,N, and /3 — 1, . . . , M. 
Let the mean of each chain be defined by 

1 ^ 

P'p = J^T.l^o.p. (43) 

Q = l 



and the mean of the distribution by 



N M 

p'=j^j:j:p^^- (44) 

a=l p=l 



We then define the variance between chains and the variance within chains as 

M 



= J^^Y.^P'^-P'f^ (45) 

13=1 



= wrhi^mi^^^^^-p'^f- (46) 

^ ' 13=1 a = l 

Lastly, we can define an estimate of the variance of p' using 

TV - 1 1 

r = __Ty« + _B^ (47) 

The estimate V"^ is an unbiased estimator of the variance if the chains are stationary, but overestimates the variance 
if the starting points of the chains are dispersed. On the other hand, typically underestimates the variance, since 
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the finite chains would not have had enough steps to range over the entire distribution. Both and approach 
the true variance in the hmit N —^ oo, but from opposite sides. If we now form the ratio 

- ^' (48) 

this ratio should approach 1 from above as the chains grow long, i.e., as the chains more faithfully reproduces the 
target distribution, and small values of i?' indicate that the chains have mixed well and converged. We take the 
condition i?* < 1.1 for all parameters as our convergence criterion. 



IV. CURRENT BOUNDS FROM CMB DATA 

We first present the constraints on the oscillation parameters using the WMAPext data set 0|, which consists of 
the WMAP 1-year data set IsS] supplemented by the Arcminute Bolometer Array Receiver (ACBAR) li^ and Cosmic 
Background I mag er (CBI) [43 data sets to extend the multipole range. Although this analysis is similar to those 
presented in j22l I2M l24| , our parametrization of the primordial power spectrum removes some of the degeneracies 
present in the original parametrization. 

In our analysis, we focus on the constraints that can be placed on the parameters defining the oscillations in the 
presence of cosmological parameters. We vary the following four cosmological parameters {ilbh^ ,^lci)Mh^ , Hq, Zrc), 
as well as the six parameters {As, Us, A, e, uj, q) that define the primordial power spectrum, assuming flatness of the 
universe. Here, As parametrizes the normalization of the primordial power spectrum according to 

5l - 10-i°A„ (49) 

and matches the normalization convention used in CAMB. 

The likelihood contours are calculated using CosmoMC, with the power spectra computed from a version of CAMB 
modified to incorporate the oscillatory power spectrum. We assume a flat prior for all parameters, restricting the 
ranges of uj and e to be w < 20 and e < 0.2, respectively. 

Our results are based on eight chains started at random points in the 10-dimensional parameter space. We flnd 
that the chains have converged after 100000 points, using the convergence criterion presented in Sec. IIII Bl 

The limits on the oscillation parameters are given in Table ^ with the corresponding likelihood contours shown 
in Fig. 121 Not surprisingly, the current data set does not give a strong constraint. At 2a, the amplitude of the 
oscillations are limited to A < 0.77, indicating consistency with a scale-invariant power spectrum. Assuming \X\ = 1, 
this translates to a restriction of the cutoff scale to A > l.SiJinfl, and does not lead to any conclusions if we allow 
|X|<1. 

The lack of a strong constraint on A is largely due to a large degeneracy between three parameters for a nonoscillatory 
fiducial model. Specifically, the three subspaces with A ^ 0, w ^ 0, and e ^ all correspond to lack of oscillations, 
so that a data set consistent with a lack of oscillations will have flat likelihoods along these planes. This behavior can 
be seen in the likelihood contours of Fig. |21 where we have marginalized the likelihood over all parameters except A 
and Lo. As can be seen, the contours curve along the A = and w = lines, indicating that both A = and lo = {) are 
supported, so that strong constraints cannot be placed on either of the two parameters. 

Previous constraints using the current data set come from Refs. |2^l23l |. Reference [22] flnds that the current data 
set is consistent with a scale-invariant power spectrum, only excluding the the region with A > 0.1 and 0.1 < e < 1. 
Since our results indicate that, for the A-e plane, the exclusion region is given by A > 0.5, e > 0.05 to 2cr, the two 
analyses are roughly in agreement. 

In contrast, Ref. !23| flnds that a model with (A = 0.29, uj = 288, e = 0.02) has an effective that is lower than the 
standard model without oscillations, with a greater than 3cr signiflcance. We deliberately avoid the high-frequency 
regime, since gravitational lensing tends to suppress the oscillatory features for models with lo > 30. We therefore 
cannot make any statements regarding the consistency of the results |47|. 

Lastly, we comment on the relation between the cosmological parameters and the oscillation parameters. As 
expected, we find that the cosmological parameters are essentially independent from the oscillation parameters. The 
cosmological parameters (e.g., fif,/i^) infiuence the structure of the acoustic peaks, which have roughly equal spacings 
in k. On the other hand, the oscillations we consider are in log(fc/fc,), and intrinsically vary separately from the 
acoustic peaks. As an illustration, we show the likelihood contours for flcBuh^ and e in Fig. ^ Thus, the additional 
parameters do not bias the measurements of the cosmological parameters. 
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V. 



THEORETICAL LIMIT 



We are ultimately interested in the intrinsic limitations of the CMB as a probe of primordial physics. To investigate 
this, we now consider an experiment that is noiseless and measures the multipole moments out to some maximum 
multipole /max: so that the uncertainties in the measurement are due to cosmic variance alone. With this idealized 
experiment, we use the transfer functions from the hierarchy code and explore the six-dimensional parameter space 
defining the primordial power spectrum with high precision. Due to the computational expense required, we choose 
to keep the cosmological parameters fixed. While such a situation may be unrealistic, the resulting limit on our 
ability to measure the parameters of the primordial power spectrum is robust, in that any consideration of realistic 
effects can only make the uncertainties larger. Furthermore, since the results from the previous section indicate that 
the cosmological parameters are independent of the oscillatory parameters to Wmax — 20, fixing the cosmological 
parameters should not bias our results. Thus, the limitations on parameters we quote should be considered as 
theoretical limits on the detectability of cutoff physics, arising purely from the statistical properties of the microwave 
background and the processing due to the radiation transfer functions. 

The maximum multipole /max of the idealized experiment is set by considering the additional complications due to 
gravitational lensing and other secondary effects. The power from gravitational lensing becomes comparable to that 
from the primary anisotropics at I 3000, and presents a fundamental limit to the range of multipoles we can use |33j . 
However, other secondary effects, such as the thermal and kinetic Sunyaev-Z'eldovich effect become dominant by 
I ~ 3000. Furthermore, deep field observations by CBI suggest an excess of power in the range I « 2000 — 3500, which 
may be due to the aforementioned secondary effects • We therefore choose a conservative value of /max — 2000 to 
avoid the regime where secondary effects may be dominant. 

For this ideal case, we simulate a data set by generating a realization of CMB power spectra that corresponds to a 
primordial power spectrum with = 5.07 x 10^^, = 0.99, and without the oscillatory features. For convenience, 
we reparametrize the normalization to 



so that Af = 1 corresponds to the fiducial model. This is related to the normalization As hy Af = 3.89 x 10~^As, 
and is chosen for computational convenience. The cosmological parameters are fixed to be {i7f,/i^ — 0.024, ilmh^ = 
0.1399, Ude — 0.73, t = 0.17, T/S — 0, wde ~ —1} for both the simulated data set and the model power spectra. From 
a converged run consisting of three chains with 230000 points per chain, we find the maximum likelihood point and 
uncertainties on the six parameters {Af^ Us, A, e, uj, a} that specify the primordial power spectrum, where we impose 
the priors e < 0.2 and cj < 30 to avoid the regime where gravitational lensing becomes important. Gravitational 
waves are assumed to be absent. 

Since we have taken the oscillations to be absent in our simulated data set, the resulting constraints on the oscillation 
parameters can be regarded as the extent to which the oscillatory parameters can deviate from the fiducial model 
before being statistically different from the standard case. For instance, a constraint of A < 0.1 at Icr would imply 
that models with A < 0.1 are statistically indistinguishable from the standard model with A = 0. Such a bound can 
be thought of as the fundamental lower limit achievable by a CMB experiment. 

We present the limits in Table ^ and the likelihood contours in Fig. [SJ showing that a cosmic variance limited 
experiment to / k, 2000 can limit the oscillations to an amplitude of less than 0.5% and frequencies of a; > 6.5, 
although other parameters are less constrained. As noted in Ref. '22], the likelihood function has many local maxima 
and irregularities. Our parametrization seems to reduce the irregularities in the amplitude parameter A; however, 
features are evident in the two-dimensional likelihood contours of the other parameters. While we believe our chains 
have converged, the presence of such features in the contour plots suggest that robust predictions for parameter 
uncertainties may require calculation of the likelihood function on much finer intervals. 

In particular, a plot of the likelihood function over e and lo, with fixed amplitude and phase, reveals that the local 
maxima are distributed along lines of constant w/e, and become sparse between w = 4 and w = 5 (Fig. The 
distribution along constant u)/e can be understood as follows. Given that the oscillation for a = is given by 



shifting oj by 2n7re/ (k/k^Y for some integer n does not change the value of the cosine at wave number k. Furthermore, 
{k/k^Y is a slowly varying function of e, so that the value of the cosine is roughly constant for uj / 1 in integer multiples 
of 27: / {k / k^:Y . If we now make a simplifying assumption that the likelihood is dominated by contributions from a 
small range in k around some value kd, then the likelihood function should exhibit this multiplicity as structures along 



5} = (5.07 X IQ-'^fAf, 



(50) 




(51) 
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The low likelihood in the range 4 < w < 5 can be explained by noting the frequency of oscillatory peaks and troughs 
in each decade. In each decade, there are approximately wln(10)/27r oscillations, so that in the range 4 < w < 5, this 
corresponds to roughly three extrema in the oscillatory features (see Fig. Cj). The extrema in the oscillations tend to 
coincide with the CMB acoustic peaks, and the overall effect of the primordial oscillations is to change the heights of 
the acoustic peaks. As is well known, the ratio of heights of the CMB acoustic peaks can be constrained tightly, and 
as a result, the likelihoods for 4 < < 5 tend to be lower. Because our MCMC sampler moves the chain based on 
the ratio of likelihoods, the lack of high likelihood regions between w = 4 and lu = 5 implies that the chains do not 
tend to cross into this region, possibly explaining the fact that our chains do not explore the region below to = 5. 

Given the presence of such irregularities, obtaining reliable constraints on the oscillation parameters is difficult. 
We believe that our reparametrization of the oscillation parameters into physically distinct parameters reduces the 
irregularities in the amplitude parameter, making the results concerning A robust. However, to obtain robust limits 
on the other parameters, the presence of multiple maxima needs to be taken into account. 

In MCMC analysis, the step sizes are typically taken to be similar to the characteristic size of the (hopefully global) 
maximum of the likelihood. However, if the local maxima have similar sizes, the chains can get stuck in the local 
maxima for a long time and lead to false conclusions about the location and constraints from the chains. In exploring 
these models, it is therefore important to monitor convergence using multiple chains, since it is unlikely that chains 
started from different, randomly chosen starting points would all fall into the same false local maximum. Furthermore, 
nonstandard algorithms (e.g., high-temperature chains and other modified samplers ^3S|) might be necessary to 
improve mixing of the chains, allowing exploration of the irregular features in a reasonable amount of computation 
time. 



VI. DISCUSSION 

We have performed a full likelihood analysis of the detectability of oscillatory primordial power spectra in current 
and future CMB data, motivated by a phenomenological description of possible effects on the primordial power 
spectrum due to the presence of a high-energy cutoff (2l| . Not surprisingly, only a weak constraint can be placed on 
the parameters that define the oscillations, the current data set being consistent with a scale-invariant power spectrum 
without oscillations. 

The analysis of the fundamental limitation from a cosmic variance limited experiment gives us cause for some 
cautious optimism, suggesting that the amplitude of the oscillations can be limited to less than 0.5%. Under the 
assumption that the vacuum choice parameter \X\ « 1, this translates to a possible limit on the cutoff energy scale of 
A > 200-ffinfi. Stated differently, if A < 200-ffinfl, the resulting features in the primordial power spectrum is observable 
in principle, allowing such a model to be distinguished from the null model. Note that this is currently the only 
limit on the cutoff scale using the effects on the vacuum choice. As noted above, however, the limits on the other 
parameters may not be robust due to the fact that the likelihood function exhibits multiple local maxima. Future 
analyses may need to take this into account by utilizing MCMC techniques that facilitate jumps between maxima, 
such as modifications to the sampling method. 

Is it possible to do better? At present, the CMB remains the best tool for exploring the shape of the primordial 
power spectrum, limited fundamentally by gravitational lensing due to the intervening matter 15]. Other probes of 
cosmological fluctuations, such as weak lensing surveys, can help extend the wave number range that can be probed, 
but exploring the effect of uncertainties from such surveys is beyond the scope of this work. 

Given the complex structure of the likelihood function, robust results from high-precision surveys will require the 
use of modified MCMC algorithms in order to map out the multiple features in the likelihood. With that caveat, 
future surveys, such as the Planck and CMBPol missions, should be able to limit the shape of the primordial power 
spectrum and provide an insight into the nature of the high energy cutoff. 
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Tables 



Parameter la 2a 
A < 0.20 < 0.77 

€ - - 

< 10.1 < 17.8 

a - - 

TABLE I: Constraints on A^(fe) parameters from the WMAPext data set, indicating the la and 2a uncertainties on parameters. 
Dashes indicate that no meaningful constraints were possible. 



Parameter 


mean 


la 


2a 


A 


0.002 


A < 0.003 


A < 0.005 


e 


0.10 


0.07 < e < 0.14 


0.002 < e < 0.19 




18.5 


16.1 <U3 < 22.3 


6.8 < w < 28.6 


a 


3.6 


2.6 < a < 4.7 


0.7 < a < 6.0 



TABLE II: Constraints on A^(/c) parameters for a cosmic variance limited experiment, assuming a fiducial model without 
oscillations. In particular, the amplitude of oscillations can be constrained to be below 0.5%. 



Figures 




FIG. 1: Fractional deviation of CAMB computation from the true power spectrum, for (a) A = 0.01, e = 0.2, co = 20, a = 0, 
and (b) A = 0.01, e = 0.2, a; = 30, q = 0. At lo = 20, the deviations are at around 0.2%, while at u = 30, the deviations 
approach 0.4%. 
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FIG. 2: WMAP likelihood contours for parameters describing the primordial power spectrum. The diagonals show the 
marginalized likelihoods as functions of the individual parameters. The off-diagonals show the two-dimensional likelihoods 
after marginalizing over the other parameters. Solid lines denote the la contours, and the gray lines the 2a contours. 




FIG. 3: Marginalized likelihood plot for the parameters A and with the solid curve denoting the la contour and the gray 
curve the 2a contour. The contours curve around the A = and uj = lines and exhibit the degeneracy between A = and 
w = models. 
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FIG. 5: la and 2a likelihood contours from a cosmic variance limited experiment, with a simulated data set containing no 
oscillations. The gray curved denote 2a contours. The diagonals show the marginalized likelihoods as functions of the individual 
parameters. 
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FIG. 6: The likelihood function, computed on a grid of u and e, holding the other variables fixed at A = 0.002, a = 0. The 
curves are iso-likelihood contours, the levels evenly spaced between minimum and maximum, with the darker spots having 
higher likelihood. Near w = 5, the local maxima seem to be lower than elsewhere, possibly accounting for the fact that the 
MCMC rules out w < 6. 




FIG. 7: Ratio of the oscillatory Cis to a fiducial Ci without oscillations, plotted for three values of e. Here, ui = 4.4, and the 
phase is fixed to be zero. The solid curve denotes Cfs with e = 0.08, the dotted curve denotes e = 0.1, and the dashed curve 
denotes e = 0.12. Note the overall coincidence of the extremum near I = 200, where the first acoustic peak lies. 



